Dilution Effects in Two-dimensional Quantum Orbital System 
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We study dilution effects in a Mott insulating state with quantum orbital degree of freedom, 
termed the two-dimensional orbital compass model. This is a quantum and two-dimensional version 
of the orbital model where the interactions along different bond directions cause frustration between 
different orbital configurations. A long-range correlation of a kind of orbital at each row or column, 
termed the directional order, is studied by means of the quantum Monte-Carlo method. It is shown 
that decrease of the ordering temperature due to dilution is much stronger than that in spin models. 
Quantum effect enhances the effective dimensionality in the system and makes the directional order 
robust against dilution. We discuss an essential mechanism of the dilute orbital systems. 

PACS numbers: 71.10.-w, 71.23.-k, 75.30.-m 
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Recently, there is much effort to reveal impurity effects 
on correlated electron systems. One of the famous exam- 
ples is the non-magnetic impurity doping into the low- 
dimensional quantum-magnets, such as the spin-Peierls 
systems and even-leg spin ladders, where the antiferro- 
magnetic (AF) long-range order (LRO) is induced by 
breaking the quantum spin-singlet correlation. 

Randomly dilute correlated electron systems with or- 
bital degree of freedom is another class of disordered 
systems. Since this degree of freedom corresponds to 
anisotropy of the electronic clouds, the orbital interaction 
in a solid is essentially directional. That is to say, Hamil- 
tonian explicitly depends on the bond direction unlike 
spin models 0, Q • There are competition and frustration 
between the different-bond interactions which reduce the 
effective dimensionaliry in the system. Dilution effects 
in such a unique system have been studied experimen- 
tally in KCui-^Zn^Fa One of the doubly degenerate 
e g orbitals in a Cu ion is occupied by a hole and ex- 
hibits LRO in KCuF 3 . With substitution of Cu by Zn, 
the orbital ordering temperature Too rapidly decreases 
and disappears around x — 0.5. This decrease of Too 
is much stronger than that in dilute magnets and in the 
percolation theory. 

This experimental observation was well explained by 
the calculation based on the orbital Hamiltonian [jj: 



H = 2J 
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with a factor e ri taking one for Cu and zero for Zn, and 
the bond direction I. Here T ri is the pseudo-spin (PS) 
operator at site i, and rh l is a unit vector in the T x — T z 
plane. When we take m l = [sin (^tt) ,cos (^-tt)] with 
(n x ,n yi n z ) = (2,1,0), this model represents the inter- 
action between the e g orbitals. Essence of the orbital 
dilution is frustration between the different-bond inter- 
actions and a macroscopic orbital degeneracy in the clas- 
sical ground state due to the interactions. 

In this paper, we study dilution effects in the two- 
dimensional and quantum version of Eq. |TjQ| with m x = 



[1,0] and rh z = [0,1], termed the quantum orbital- 
compass model (OCM) 0, Q. This is a minimal or- 
bital model representing the directional and frustrated 
interaction between orbitals. One of the major different 
points from the e s -orbital model is that the Hamiltonian 
is invariant under a certain local transformations of the 
PS operator. Due to this symmetry, degeneracy of the 
orbital state is more serious than that in the e ff -orbital 
model, and the conventional site-diagonal order is not ex- 
pected. Instead, the long-range correlation for one of the 
PS components appears along the vertical or horizon- 
tal direction in a square lattice, termed the directional 
order (DO). Dilution effects in DO are investigated by 
means of the quantum Monte-Carlo (QMC) method in 
addition to some rigorous analytical considerations. We, 
in particular, focus on the quantum effects in the orbital 
dilution. It is shown that quantum fluctuation from the 
classical DO state enhances the effective dimensionality 
in the system and makes DO robust against dilution. 

The explicit form of OCM is given in a two-dimensional 
square lattice as 



H = 2J 
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Since the sign of the interaction J is gauged away by 
rotating T ri with respect to the y axis, we take it neg- 
ative. This Hamiltonian without dilution (e ri = 1 for 
Vi) is invariant under the following two-symmetry op- 
erations 0, II] ; (i) The global four-fold symmetry: the 
pseudo-spins at all sites and a crystal lattice are rotated 
by 7r/2, simultaneously, with respect to the y axis, (ii) 
The local symmetry at each column and row: the z (x)- 
component of all PSs at each column (row) along the 
z (x) axis are flipped, i.e. T^) z — > —Trj; x J z for each 
r x ( r z)- In addition, in the classical ferromagnetic- like 
ground state, there is the following continuous symmetry: 
(iii) the global PS direction is arbitrary in the T x — T z 
plane. 

Let us explain in more detail the symmetry (ii) which 
is unique in OCM. This symmetry operation at a row 
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is done by the operator P r 
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the Pauli matrices a 
-l 



with 

as Pr z la %^,r z Pr z = — (X r x ,r z 

Pr} a r x .r,Pr z = o* r . O n an equal footing, the sym- 
metry operation at a column r x is done by the operator 
Q rx — Yl r of r ■ The operators P rz and Q ri commute 
with Tt, and P rz s commute with each other as well as 
Qr x s, but [Qr x ,Pr x ] 7^ 0. The energy eigen- values are 
characterized by the eigen- value set of P rz s (pi, ■ ■ ■ ,Pl) 
or that of Qr x s (qi, ■ ■ ■ , q£) in a L x L(= N) square lat- 
tice. That is, at least, there are the L conserved quan- 
tities. Since P% ( = 1 and pt = P n , pi = ±1. It 
is easily shown that the conventional site-diagonal or- 
der parameter N _1 5Z r (T r ) is not invariant under this 
local-symmetry operation. According to the generalized 
Elitzur's theorem, it is shown that the site-diagonal or- 
bital LRO does not occur 0. Instead, a certain kind 
of LRO, termed the directional order, appears. This 
is a long-range correlation of T z (T x ) at each column 
(row), and corresponds to a breaking of the global four- 
fold symmetry (i) explained above. This nematic-likc 
order has been studied by means of the correlation func- 
tion Cy* = hm rs ._Kx, (T^ T^ r ^ ) at zero temperature la 
and q = iV" 1 J2r { T f ~ T^J in the classical OCM 
We introduce, as an order parameter of DO, a more con- 
venient form D defined by 
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This parameter D indicates the correlation of T z along 
the z axis or that of T x along x. Unlike q, this is not a 
constant for the quantum PS operators. Since D com- 
mutes with P Tz and Qr x , being different from T r , it is 
possible for (D) to be finite. 

To examine OCM numerically with and without va- 
cancies, we adopt QMC method in a finite-size cluster. 
There is no negative-sign problem in OCM. The QMC 
calculations have been performed in a square lattice of 
L x L sites (L=14~18) with the periodic-boundary con- 
dition. We chose the Trotter number n to be 12~22, 
and extrapolate the calculated results to n — oo. At low 
temperatures of T < 0.04 J, wc perform the simulation 
with n =60^72. We take, at maximum, 10,000MC steps 
for measurement after 3,000MC steps for thermalization. 
Physical quantities are averaged over 20~60MC samples 
at each parameter set. We have checked by MC calcula- 
tions with 10 6 MC steps, which is much longer than the 
integrated autocorrelation time determined by the bin- 
ning analyses, that the order parameter values in the 
present calculation are reliable within numerical errors. 
The ordering temperature (Tbo) of DO is estimated from 
the finite-size scaling analyses of the Binder cumulant for 

D defined by Q = 1- 0,E3- At T = T DO , a value 

of Q should be size independent. 

First we show the numerical results in OCM without 
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FIG. 1: (color online), (a) Temperature dependence of the 
directional order parameter D(T) at x = for various system 
sizes L. The inset shows D(T) at low temperatures obtained 
in the calculation with the large Trotter numbers (filled cir- 
cles). Open circles are the same data with those in the main 
panel, (b) D(T) for various impurity concentrations. 



vacancies. Temperature dependence of y/ (D 2 )(= D(T)) 
for several L is presented in Fig. 1(a). Around T/J = 
0.17, D(T) grows with decreasing T, and is saturated to 
about 0.13 below T/J = 0.12. With increasing L, this 
dependence becomes steep, but the saturated value of 
D(T) at low temperatures does not change much. Be- 
low T/J = 0.07, D(T) seems to decrease. We have 
performed careful numerical calculations with the larger 
Trotter numbers n =60^72. As shown in the inset of 
Fig. 1(a), ~D(T) takes about 0_T25 down to T/J = 0.015. 
Therefore, the reduction of D(T) in low temperatures 
is an artifact due to the small Trotter number. Tem- 
perature dependence of the Binder cumulant Q is pre- 
sented in Fig. 2. The Q — T curves for the several sys- 
tem size L cross around T/J = 0.15 corresponding to 
the temperature where D(T) grows. We have also calcu- 
lated the susceptibility xd of D and the specific heat C. 
Both of them show peak structures around T/J = 0.15, 
which become steep with increasing L. We conclude in 
the present calculation that DO is realized in the quan- 
tum OCM and T DO =^_0.150 ± 0.003. It is noted that 
the saturated value of D at low temperatures is about a 
half of the classical value of 1/4. This reduction is much 
stronger than that in the quantum AF Heisenberg model 
with 5 = 1/2 in a square lattice where the saturated 
staggered moment is about 60% of 1/2 ^lj- We also cal- 
culate the site-diagonal order parameter M = yj (M 2 ) 
with M = N^ 1 J2 r Tr- The calculated M is almost tem- 
perature independent and its magnitude decreases with 
increasing L. This is, the site-diagonal order does not 
tend to occur because of the local symmetry (ii), men- 
tioned previously. This situation is in contrast to the 
e s -orbital model where the site-diagonal order is realized 
by the order-by-disorder mechanism. 

Now we present the dilution effects in OCM. D(T) for 



3 



0.6 

a 

0.55 

0.5 

0. 



• i 



x=0 



L • 

• 14 ♦ 

* 16 ■■. 

■ 18 *' 

J i l_l 



L 

• 14 

* 16 
■ 18 

f 

■ * 



14 



0.15 0.16 
T/J 



FIG. 2: (color online) . Temperature dependence of the Binder 
cumulant Q at x — for various system sizes L. Inset shows 
temperature dependence of Q at x = 0.06. 



several x are shown in Fig. 1(b). D(T) gradually in- 
creases with decreasing T and is saturated, for example, 
to about 0.1 at x = 0.06. The temperature where D(T) 
rises up is reduced with doping. The saturated value 
of D(T) becomes small by doping in spite of a factor 
(1 — x)^ 1 in Eq. (J3J. This phenomenon was also ob- 
served in the dilute e ff -orbital model j|| and will be dis- 
cussed in more detail later. The analyses for Q work well 
at x = 0.06, as shown in the inset of Fig. 2, where Tdo / J 
is identified to be 0.113 ± 0.005. We have checked that, 
at x = 0.25, the Q — T curves for several L do not cross 
with each other down to T = 0.01J. 

The main results in the present work, the x dependence 
of Too in the quantum OCM, are shown in Fig. 3. For 
comparison, we plot the ordering temperature Tcl of DO 
in the classical OCM and the Curie temperature Ti s ; ng 
in the 5 = 1/2 Ising model in a square lattice. Tdo 
decreases monotonically with x and seems to take zero 
around x = 0.15(= x c ). It is worth to note that x c is 
much smaller than the percolation threshold 0.41 in a 
square lattice |4J . The decrease of Tdo is much stronger 
than that in Ti s i ng , but is weaker than Tcl- In other 
words, quantum effect makes DO robust against dilution. 

To understand the unusual dilution effects in the quan- 
tum OCM, we present a rigorous treatment of the di- 
lute OCM in two-coupled chains along the z axis, that 
is, a two-leg ladder OCM. As mentioned previously, the 
energy eigen-states are classified by the eigen-value set 
(l>i) •'• >.£>-£/) °f the operator P Tz — a* T a* r where 
r R( r L) indicates the x coordinates of the right (left) chain 
in the ladder. We consider one of the two-fold ground 
state with pi = ■ ■ ■ = pi, = 1 being degenerate with the 
state characterized by pi = ■ ■ ■ = pl = — 1. In a row r z 



of the ladder, the eigen states are \T 1 
and | ||), or their linear combinations |±) 
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FIG. 3: (color online). Impurity concentration x dependence 
of the ordering temperature. Closed circles are for Tbo in 
the quantum OCM, and closed squares are for Tcl in the 
classical OCM. For comparison, the ordering temperature in 
the Ising model (triangles) are also plotted. Bold arrow in- 
dicates the percolation threshold x c in a square lattice. The 
inset shows the ordering temperature Tc in the quantum XY 
model (circles) and that in the classical one (squares) in a 
three-dimensional cubic lattice. System size is chosen to be 
TV = 23 3 at maximum. The continuous imaginary-time algo- 
rithm is utilized. 
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\±) = 



Therefore, 



within the subspace of pi = • • • = pl = 1, we intro- 
duce new pseuso-spin operators S rz and show correspon- 
dences (T^T^J - and 2^ (rjl)r , - S* . The 

state |+) (|— )) is the up- (down-) eigen state of the op- 
erator S^. As a result, OCM in this subspace is mapped 
onto the transverse-Ising model in a single chain |13| : 
H e ff = ~^JJ2r z Sr ss Sr ss +i~JYlrz &r z - Now one vacancy 
is introduced at a site (rL,r"). Since, in the row r°, we 
have P r o — <Jr Rr o, the energy eigen states in this row, 
|— t) and |— |), belong to the different subspace, and 
T rRT o is replaced by a C-number in each subspace. Thus, 
OCM in the two-leg ladder with one vacancy is denoted 
by the following model: 

/ 

± JiS^ + S^), (4) 

where ^2' r indicates a sum of r z except for and r® — l. 
Signs + and — in the third term correspond to the states 
|— T) and |— |) at the site {rR,r®), respectively. This 
model is the two transverse-Ising chains where the exter- 
nal field ± J along the S x axis is applied at edges of the 
chains, i.e. the r° + 1 and r° z — 1 sites . In the original 
language, this field acts on T* n{rL ^ r o +1 and T* R ^ r ^ r0 _ v 
and enhances the correlation (T*.T*.) around the va- 
cancy. The above analyses accord with an intuitive pic- 
ture for roles of vacancy in OCM; due to the directional 
interaction in the model, orbitals near the vacant site ad- 
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FIG. 4: (color online) . Contour map of local static correlation 
function C zz (ri) in a 20 x 20 site cluster . Four vacancies 
Xi(2) and Zi(2) are represented by black circles. Temperature 
is chosen to be T/T DO (x = 0) =1.27. 



just their directions to gain the exchange energy. In other 
words, vacancy tends to fix the orbital direction around 
it. 

Being based on the above picture, we discuss the dif- 
ference between Tbo and Ti s i ng in Fig. 3. Consider a 
situation in a two-dimensional square lattice where two 
vacancies are located closely in a same row, e.g. at sites 
(r x — 1 , r z ) and (r x + l,r z ). This probability is of the order 
of x. Since PS between the vacancies has the Ising-type 
interaction T* rx (T* xrz+1 + T z ^ r _ 1 ), the correlation of 
T z is enhanced around the vacancies, as suggested above. 
On an equal footing, the correlation of T x is enhanced 
around the vacancies located closely in a same column. 
It is expected that a certain number of domains/patches, 
where T z or T x correlation is enhanced, coexist in a sys- 
tem. Competition between such kinds of domains reduces 
the order parameter D(T) and the ordering temperature 
Tdo- To confirm this picture, we calculate the local- 
static PS correlation function around the impurity de- 
fined by C7"( r< ) = £ . (T r \T*.) (Fig. 4). Temper- 

ature is chosen to be T/Tdo{% = 0) = 1.27. Two pairs of 
vacancies introduced in a 20 x 20 site cluster are denoted 
by -X"i(2) and ^1(2) in the figure. It is clearly shown that 
C zz (vi) is enhanced between the impurities Xi and X2, 
and almost vanishes between Zi and Z2. These results 
are not seen in dilute spin models. 

Next we focus on quantum effects in dilution, i.e. dis- 
crepancy between Tdo and Tcl- In the classical DO 
state at T = 0, the system is considered to be decou- 
pled into the independent Ising chains along the 1 or z 
axis, although DO is the two-dimensional order where the 
four-fold symmetry is broken. The rapid decrease of Tex 
by dilution reflects the quasi-one dimensional nature of 
DO, and at finite temperatures, weak two-dimensionality 
is recovered due to the thermal fluctuation. The quan- 
tum fluctuation from the classical DO state brings about 
a coupling between the independent chains and increases 



effective dimensionality of DO even at low temperatures. 
This is a kind of resonant states between the independent 
chains and DO becomes robust against dilution. This re- 
markable quantum effect is unique in the orbital system; 
we compare the x dependence of the Curie temperatures 
Tq in the quantum and classical XY models in a cubic 
lattice (the inset of Fig. 3). A magnitude of spin is cho- 
sen to be 1/2 in the quantum model. It is seen that 
difference between the two are much smaller than that in 
OCM, and both Tc's in the XY model seem to disappear 
at the percolation threshold (x c =0.69). 

Before summarizing the results, we briefly discuss im- 
plications of the dilute e g -orbital model. In our previous 
studies |5j , we have suggested that a remarkable decrease 
of Too is related to the tilting of PS around impurities. 
Roles of this PS tilting become clearer by the present 
analyses of OCM. The orbital interaction in the e g -orbital 
model is essentially of the Ising-type along one direction 
in a cubic lattice, and there is frustration between them, 
as well as OCM. Introduced vacancies break the local 
symmetry and release the frustration. According to va- 
cancy configurations, a specific type of the orbital corre- 
lation is stabilized around vacancies. For example, when 
two vacancies are located closely at a same row along 
the z axis, the orbital correlation of d^ x 2_ r 2 / d y 2_ z 2 or 
d 3y 2_ r 2 /d z 2_ x 2 are stabilized around them. Competition 
between these domains/patches with different PS corre- 
lations reduce the order parameter and may cause the 
glass-like phase. 

In summary, we study the dilution effects in a quantum 
and two-dimensional version of the orbital model, i.e. the 
orbital compass model, where the vertical or horizontal 
correlation of a PS component, termed the directional 
order, is expected from the local symmetry. Reduction 
of the ordering temperature of DO is much stronger than 
that in the spin model. Quantum fluctuation increases 
the effective dimensionality of the system and makes DO 
robust against dilution. Competition between domains 
with different orbital configurations is essence of the di- 
lute orbital system. 
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